Astronomy & Astrophysics manuscript no. alphacentauri 


© ESO 2008 


June 17, 2008 





The Alpha Centauri binary system * 
Atmospheric parameters and element abundances 

G. F. Porto de Mello, W. Lyra **, and G. R. Keller'*** 

Observatorio do Valongo, Universidade Federal do Rio de Janeiro, Ladeira do Pedro Antonio, 43 
20080-090 Rio de Janeiro, RJ, Brazil 

Received, accepted 

ABSTRACT 

Context. The a Centauri binary system, owing to its duplicity, proximity and brightness, and its components' likeness to the Sun, is a funda- 
mental calibrating object for the theory of stellar structure and evolution and the determination of stellar atmospheric parameters. This role, 
however, is hindered by a considerable disagreement in the published analyses of its atmospheric parameters and abundances. 
Aims. We report a new spectroscopic analysis of both components of the a Centauri binary system, compare published analyses of the system, 
and attempt to quantify the discrepancies still extant in the determinations of the atmospheric parameters and abundances of these stars. 
Methods. The analysis is differential with respect to the Sun, based on spectra with R = 35 000 and signal-to-noise ratio > 1 000, and employed 
spectroscopic and photometric methods to obtain as many independent T e g determinations as possible. The atmospheric parameters are also 
checked for consistency against the results of the dynamical analysis and the positions of the components in a theoretical HR diagram. 
Results. The spectroscopic atmospheric parameters of the system are found to be T e s = (5847 ± 27) K, [Fe/H] = +0.24 ± 0.03, log g = 4.34 
± 0.12 and £ t = 1.46 ± 0.03 kms" 1 , for a Cen A, and T eff = (5316 ± 28) K, [Fe/H] = +0.25 ± 0.04, logg = 4.44 ± 0.15 and f t = 1.28 ± 0.15 
kms -1 for aCenB. The parameters were derived from the simultaneous excitation & ionization equilibria of Fel and Fell lines. T e gs were 
also obtained by fitting theoretical profiles to the Ha line and from photometric calibrations. 

Conclusions. Good agreement was reached between the three criteria for a Cen A. For a Cen B the spectroscopic T e g is ~ 140 K higher than 
the other two determinations. We discuss possible origins of this inconsistency, concluding that the presence of NLTE effects is a probable 
candidate, but we note that there is as yet no consensus on the existence and cause of an offset between the spectroscopic and photometric 
T e ff scales of cool dwarfs. The spectroscopic surface gravities also agree with those derived from directly measured masses and radii. An 
< £j ■ average of three independent T eB criteria leads to T eff (A) = (5824 ± 26) K and T eff (B) = (5223 ± 62) K. The abundances of Na, Mg, Si, Mn, 
k> Co and Ni and, possibly, Cu are significantly enriched in the system, which also seems to be deficient in Y and Ba. This abundance pattern can 
» I ■ be deemed normal in the context of recent data on metal-rich stars. The position of a Cen A in an up-to-date theoretical evolutionary diagrams 

' yields a good match of the evolutionary mass and age (in the 4.5 to 5.3 Gyr range) of with those from the dynamical solution and seismology, 

_ _ i 

but only marginal agreement for aCenB, taking into account its more uncertain T c g. 
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1. Introduction 

The a Centauri binary system, composed of two solar-type 
stars (HD 128620 and 128621), is one of the brightest in the sky 
and figures as our second closest galactic neighbor, 1.34 parsec 
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away. The star closest to the Sun is the M5.5 dwarf Proxima 
Centauri (Gliese & Jahreiss 1991), ~15,000 A.U. away from 
the a Centauri binary, and its gravitational connection to the 
system is still a topic of controversy. Anosova et al. (1994) 
proposed that Proxima has a hyperbolic orbit around the in- 
ner pair, and that the three stars might form part of a more 
extended kinematical group. Wertheimer & Laughlin (2006), 
however, found the distance between Proxima and the pair as 
comparable to the Hill radius of the latter, whereby the Galactic 
potential becomes dominant over that of the inner pair and the 
system becomes unbound. These authors favor the existence 
of a physically bound triple system, suggesting that Proxima 
is presently at the apoastron of its orbit. Highly precise moni- 
toring of radial velocity variations of the system by Endl et al. 
(2001) constrains the upper limit of the mass of putative plan- 
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etary or substellar companions of the system as less than 3.5 
Jupiter masses (actually less than one Saturn mass if coplanar 
orbits are assumed). 

The proximity of the a Centauri system provides a well 
determined parallax, and its brightness allows the acquisition 
of extremely high-quality spectra. Moreover, its binary nature 
and relatively short period of 80 years enables the hypothesis- 
free accurate determination of masses (Pourbaix et al. 1999, 
2002). If we couple to these facts their being very solar-like, 
the a Centauri stars thus appear as objects of fundamental im- 
portance in the calibration of evolutionary tracks, theoretical 
isochrones and model atmospheres, hence the great interest in 
the precise determination of their atmospheric parameters, evo- 
lutionary state and chemical composition. 

The brightness of the system's components also favor the 
determination of internal structure and state of evolution by 
seismological observations. The analysis of the frequency spec- 
trum and amplitudes of both photometric and spectroscopic os- 
cillations in the outer layers of solar-type stars, driven by con- 
vection, can yield otherwise unobtainable information on inter- 
nal structure, such as the depth of the convection zone and the 
density and temperature profiles. They can also provide inde- 
pendent checks on stellar masses, ages and chemical composi- 
tion. Yildiz (2007), Eggenberger et al. (2004) and Thoul et al. 
(2003) have agreed on an age for the system between 5.6 and 
6.5 Gyr. Miglio & Montalban (2005) propose model-dependent 
ages in the 5.2 to 7. 1 Gyr interval. However, they also note that 
fixing the non-seismic observables, namely masses and radii, 
leads to an age as large as 8.9 Gyr, proposing that further seis- 
mological observations may be needed to clarify this apparent 
discrepancy between the independent observation of the os- 
cillation spectra and the directly measured masses and radii. 
The previous analysis of Guenther & Demarque (2000) favors 
a slightly higher age of ^7.6 Gyr. The masses are very well 
constrained at M A = 1.105 ± 0.007 and M B = 0.934 ± 0.006 
in solar masses (Pourbaix et al. 2002), which, along with in- 
terferometrically measured (in solar units) radii of Ra = 1 -224 
± 0.003 and R B = 0.863 ± 0.005 (Kervella et al. 2003) yield 
surface gravities (in c.g.s. units) of log gA = 4.307 ± 0.005 
and log ge = 4.538 ± 0.008, an accuracy seldom enjoyed by 
stellar spectroscopists. Altogether, these data pose very tight 
constrains on the modelling of fundamental quantities of inter- 
nal structure, such as mixing-length parameters and convection 
zone depths. 

Nevertheless, the state of our current understanding of this 
system still lags behind its importance, since published spec- 
troscopic analyses reveal considerable disagreement (Table [U 
in the determination of atmospheric parameters and chemical 
abundances, particularly for component B, though most authors 
agree that the system is significantly metal-rich with respect to 
the Sun. This fact is embarrassing, even in our modern era of 
massive surveys, since the individual study of key objects is 
necessary to quantify systematic errors which might be lurking 
inside huge databases, and cannot be reduced with large num- 
ber statistics. Indeed, considering only those analyses since the 
90s, eight performed a detailed analysis of the atmospheric pa- 
rameters and chemical composition of a Cen A (Furenlid & 
Meylan (1990), hereafter FM90; Chmielewski et al. (1992), 



Neuforge-Verheecke & Magain (1997), Allende-Prieto et al 
(2004), hereafter ABLC04; del Peloso et al. (2005a), Santos 
et al. (2005); Valenti & Fisher (2005); Doyle et al. (2005)), 
whilst five of them also performed this analysis for the cooler 
and fainter component aCenB (Chmielewski et al. 1992, 
Neuforge-Verheecke & Magain (1997), ABLC04; Valenti & 
Fisher (2005) and Santos et al. 2005). All these authors, but 
Chmielewski et al. (1992) and Santos et al. (2005), obtained 
abundances for many chemical elements other than iron. 

The analysis of FM90 for a Cen A is noteworthy in that it 
was the first to imply an abundance pattern considerably dif- 
ferent from solar, with excesses relative to Fe in Na, V, Mn, 
Co, Cu, and deficits in Zn and the heavy neutron capture el- 
ements, and also proposed a low T c s and a near solar metal- 
licity for component A, in contrast with most previously pub- 
lished figures. These authors invoked a supernova to explain 
the peculiar chemical features of the system. The next analysis 
(Chmielewski et al. 1992) sustained a high T c ff and apprecia- 
bly higher metallicity for the system, which was also obtained 
by Neuforge-Verheecke & Magain (1997). The latter authors, 
moreover, found an abundance pattern not diverging signifi- 
cantly from that of the Sun, though supporting the deficiency 
of heavy elements found by FM90. 

del Peloso et al. (2005a) and Santos et al. (2005) have both 
derived a high metallicity for the system. Doyle et al. (2005) 
added to the controversy by proposing both a low T e s and a 
metallicity not appreciably above solar for component A, as 
did FM90. Their abundance pattern is, however, solar. ABLC04 
propose for both components much lower T c ff s than previously 
found by any author. Even though their metallicity agrees rea- 
sonably with that of Chmielewski et al. (1992) and Neuforge- 
Verheecke & Magain (1997), their detailed abundance pattern 
is highly non-solar and also very different from any thus far, 
with high excesses of Mg, Si, Ca, Sc, Ti, Zn and Y. Their low 
metallicity is a result of a lower adopted T c ff , as also is the case 
for the FM90 analysis. 

Doyle et al. (2005) presented the most recent abundance 
analysis of a Cen A and obtained abundances for six elements. 
They made use of the Anstee, Barklem and O'Mara (ABO) line 
damping theory (Barklem et al. 1998 and references therein), 
which allowed them to fit accurate damping constants to the 
profile of strong lines, turning these into reliable abundances 
indicators, an approach normally avoided in abundance analy- 
ses. They found [Fe/H]=0.12 ± 0.06 for the iron abundance, 
which is in disagreement with most authors using the stan- 
dard method, although agreeing with FM90. To bring home the 
point of the existing large disagreement between the various 
published results, one needs look no further than at the last en- 
tries of Table [Tj all based on very high-quality data and state 
of the art methods. These disagreements in chemical composi- 
tion lie beyond the confidence levels usually quoted by the au- 
thors. Moreover, the dispersion of the T e s values found range 
between 300K and 400K, respectively, for component A and B. 

Pourbaix et al. (1999) finish their paper thus: "we urge 
southern spectroscopists to put a high priority on a Centauri". 
Clearly, this very important stellar system is entitled to addi- 
tional attention, fulfilling its utility as a reliable calibrator for 
theories of stellar structure and evolution, and taking full ad- 
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Table 1. A non-exhaustive review of determinations of effective temperatures T e g, iron abundances [Fe/H], logarithm of surface 
gravity \ogg and microturbulence velocities £ ( for a Cen A and B. The notation used from columns 5 to 7 stands for: a - 
Excitation equilibrium; b - Wings of Balmer lines; c - Trigonometric parallax; d - Wings of strong lines; e - Ionization equilibrium; 
/ - Photometric color indexes; direct - directly measured luminosity, mass and radius. Note that, generally, the microturbulence 
velocities have not the same zero point and cannot be directly compared. Typical errors, respectively, in T e g, log g and [Fe/H] 
are < 100 K, < 0.2 dex and < 0.1 dex, but note that these estimates usually do not include systematic uncertainties, and that not 
all authors provide error determinations, neither labor under the same definitions for them. 
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vantage of its tight observational constraints towards our under- 
standing of our second closest neighbor and the atmospheres of 
cool stars. The widely differing results of the chemical analy- 
sis also cast doubt about the place of a Centauri in the galactic 
chemical evolution scenario. The goal of the present study is 
a simultaneous analysis of the two components of the system, 
obtaining their atmospheric parameters and detailed abundance 
pattern, providing an up-to-date comparative analysis of the 
different determinations, the methods used and their results. 



This paper is organized as follows. In section 2 we describe 
the data acquisition and reduction. In section 3 we describe the 
spectroscopic derivation of the atmospheric parameters and Fe 
abundance, and compare them to other recent results from other 
techniques, discussing possible sources of discrepancies. The 
chemical composition pattern and its comparison to those of 
other authors, is outlined in section 4. Section 5 is devoted to 
the analysis of the evolutionary state of the system, and section 
6 summarizes the conclusions. 
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Fig. 1. Sample spectra of the Moon, employed as a proxy for 
the solar flux spectra, and a Centauri A and B. The nominal 
resolution is R = 35 000, and the signal-to-noise (S/N) ratio is in 
excess of 1 000. Some spectral lines measured for the spectro- 
scopic derivation of the atmospheric parameters and the abun- 
dance analysis are indicated. The top, dotted spectrum is the 
Moon's, the dark gray one corresponds to a Cen A, and the 
light gray one to a Cen B. The a Cen spectra are arbitrarily 
displaced vertically. The stronger line blocking in a Cen B is 
apparent. 

2. Observations and line measurement 

Observations were performed in 2001 with the coude spectro- 
graph, coupled to the 1.60m telescope of Observatdrio do Pico 
dos Dias (OPD, Brasopolis, Brazil), operated by Laboratorio 
Nacional de Astrofisica (LNA/CNPq). As both a Cen A and B 
are solar-type stars, the Sun is the natural choice as the standard 
star of a differential analysis. The expectation of this approach 
is that systematic errors in the measurement of line strengths, 
the representation of model atmospheres, and the possible pres- 
ence of Non-Local Thermodynamic Equilibrium (NLTE) ef- 
fects, will be eliminated or at least greatly lessened, if the stan- 
dard and the analysed object are sufficiently similar. We chose 
the moon as a sunlight surrogate to secure a solar flux spectrum. 
The slit width was adjusted to give a two-pixel resolving power 
R = 35 000. A 1 800 1/mm diffraction grating was employed in 
the first direct order, projecting onto a 24/im, 1024 pixels CCD. 
The exposure times were chosen to allow for a S/N ratio in ex- 
cess of 1 000. A decker was used to block one star of the binary 
system while exposing the other, and we ascertained that there 
was no significant contamination. The moon image, also ex- 
posed to very high S/N, was stopped orthogonally to the slit 
width to a size comparable to the seeing disks of the stars. 

Nine spectral regions were observed, centered at A A 5100, 
5245, 5342, 5411, 5528, 5691, 5825, 6128 and 6242 A, with 
spectral coverage of 90 A each. The chemical species repre- 
sented by spectral lines reasonably free from blending are Na I, 
Si I, Cai, Sci, Sen, Til, Tiii.Vi, Cri, Cm, Mm, Fei, Fen, 
Co I, Ni I, Cu I, Y II, Ba II. Additional data centered on the Ha 
spectral region, for the a Cen stars and moonlight, were se- 
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Fig. 2. A plot of Voigt-fitted W\s from the Solar Flux Atlas, 
given by Meylan et al. (1993), against our gaussian measure- 
ments, for the 84 common lines. A linear regression defines the 
correction to be applied to our measured W\ & in order to lessen 
systematic errors due to inadequate gaussian profile fitting. As 
expected for non-saturated lines, we ascertained that a linear 
relationship suffices to describe the correction. The line wings, 
fitted incompletely by gaussian profiles, start to be important 
for lines stronger than ^50mA. The best fit (solid line) and the 
one-to-one line of identity (dotted line) are shown. 

cured in 2004, using a 13.5/im, 4608 pixels CCD, integrated to 
S/N - 500 and with R = 43 000. 

Data reduction was carried out by the standard procedure 
using IRAlfl After usual bias and flat-field correction, the 
background and scattered light were subtracted and the one- 
dimensional spectra were extracted. No fringing was present in 
our spectra. The pixel-to-wavelength calibration was obtained 
from the stellar spectra themselves, by selecting isolated spec- 
tral lines in the object spectra and checking for the absence of 
blends, the main screen for blends being the Solar Flux Atlas 
(Kurucz et al. 1984) and the Utrecht spectral line compilation 
(Moore et al. 1966). Gaussian fits were applied to the cores of 
the selected lines, and pixel-A polynomial fits determined. For 
the short spectral selections individually reduced, a 2nd-order 
polynomial always sufficed, the average r.m.s. of the residuals 
being 0.005 A or better. There followed the Doppler correction 
of all spectra to a rest reference frame. 

Normalization of the continuum is a very delicate and rel- 
evant step in the analysis procedure, since the accuracy of line 
equivalent width (hereafter W\) measurements is very sensitive 
to a faulty determination of the continuum level. We selected 
continuum windows in the Solar Flux Atlas, apparently free 
from telluric or photospheric lines. We took great care in con- 
stantly comparing the spectra of the two a Cen components and 

1 Image Reduction and Analysis Facility (IRAF) is distributed by 
the National Optical Astronomical Observatories (NOAO), which 
is operated by the Association of Universities for Research in 
Astronomy (AURA), Inc., under contract to the National Science 
Foundation (NSF). 
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the Sun, to ensure that a consistent choice of continuum win- 
dows was achieved in all three objects, since the very strong- 
lined spectra of the a Cen stars caused continuum depressions 
systematically larger than in the Sun. A number of pixels was 
chosen in the selected continuum windows, followed by the 
determination of a low order polynomial fitting these points. 
The wavelength coverage of each single spectrum was in all 
cases sufficient to ensure an appropriate number of windows, 
with special attention given to the edge of the spectra. Sample 
spectra are shown in Fig. Q] As will be seen below, the errors 
of the atmospheric parameters derived directly from the spec- 
tra, and the element abundances of a Cen B, are greater than in 
a Cen A, probably due to a less trouble-free normalization of 
its strongly line-blocked spectrum, and to a better cancellation 
of uncertainties in the differential analysis. 

For the determination of element abundances, we chose 
lines of moderate intensity with profiles that indicate little or no 
blending. To avoid contamination of telluric lines we computed 
for each spectrum, using cross-correlation techniques, the dis- 
placement AA the telluric lines would show relative to their 
rest position Ao, as given by the Solar Flux Atlas. We discarded 
photospheric lines closer than 2AA from a telluric line. 

The equivalent widths were measured by fitting single or 
multiple (the latter when deblending closely spaced lines) gaus- 
sian profiles to the selected lines, using IRAF. The moderately 
high spectral resolution we chose was designed to guarantee 
that the instrumental profile dominates the observed profile, 
and therefore that purely gaussian fits would adequately rep- 
resent the observed line profiles. To test the representation of 
the solar flux spectrum by the moon, we also observed, with 
exactly the same setup, spectra of daylight and the asteroid 
Ceres. A direct comparison of the moon, daylight and asteroid 
W\s showed perfect agreement between the three sets of mea- 
surements to better than 1% even for moderately strong lines. 
This lends confidence to our determination of solar gf-values 
based on W\ measured off the moon spectra. The moonlight 
spectra was actually preferred due to its higher S/N ratio as 
compared to Ceres, for which no high-quality spectrum could 
be obtained in a reasonable exposure time. Also, daylight spec- 
tra may show W\ systematic fill-in effects by up to 4%, as a 
combination of aerosol and Rayleigh-Brillouin effects (this ef- 
fect depends on observing angle and can be eliminated or min- 
imized if care is applied, see Gray et al. 2000 for details). Even 
though no difference could be measured in our spectra, we con- 
sidered it more prudent to use the moonlight spectrum as the 
solar proxy: it should be an accurate representation of the solar 
flux spectrum in the visible. 

Even at our not-so-high resolution, lines stronger than ^50 
mA begin to develop visible Voigt wings. To account for this 
effect, we performed a linear regression of our gaussian moon 
W\s against the Was measured off the Solar Flux Atlas by 
Meylan et al. (1993). These authors fitted Voigt profiles to a 
set of lines deemed sufficiently unblended to warrant the mea- 
surement of their true W\S, and they should be a homogeneous 
and high-precision representation of the true line intensities. 
We then determined the correction necessary to convert our 
measurements to a scale compatible with the Voigt-fitted Was. 
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The result is shown in Fig. 12 where the excellent correlation, 
with very small dispersion, is seen. The correction derived is 

^Voigt = ^ Q48 ± o 013 )T^moon (1) 

The r.m.s. standard deviation of the linear regression is 
2.9 mA, regarding the Voigt W\ as essentially error-free as 
compared to our data. This regression was applied to all our 
W\ measurements. We take 2.9 mA as the la uncertainty of 
our internal W\ measurements. 

3. Atmospheric parameters and Fe abundance 

A solar gf-value for each spectral line was calculated from a 
LTE, 1-D, homogeneous and plane-parallel solar model atmo- 
sphere from the NMARCS grid, as described by Edvardsson 
et al. (1993) (see http://marcs.astro.uu.se, and also Gustafsson 
et al. 2008). The adopted parameters for the Sun were T c g= 
5780K, log g = 4.44, [Fe/H] = +0.00 and & = 1.30 kms" 1 , 
and we employed the Was measured off the moon spectra, 
corrected to the Voigt scale. The adopted solar absolute abun- 
dances are those of Grevesse & Noels (1993). In a purely dif- 
ferential analysis such as ours, the absolute abundance scale is 
inconsequential. We provide in Table [2] the details of all lines 
used. They include wavelength A, excitation potential x> the 
calculated solar log gf values, and the raw measured W\S in 
the moon's, a Cen A and aCenB spectra, prior to the correc- 
tion to the Voigt system (Fig. |2}. Hyperfine structure (HFS) 
corrections for the lines of Mg I, Sc I, Sc II, V I, Mn I, Co I, 
Cul and Ball were adopted from Steffen (1985). del Peloso 
et al. (2005b) discuss the influence of adopting different HFS 
scales on abundance analyses of Mn and Co, concluding that it 
is small, particularly for metallicities not too far from the so- 
lar one, as compared to not using any HFS data. Therefore, the 
source of the HFS corrections is not an important issue on the 
error budget of our analysis, at least for Mn and Co. The other 
elements of our analysis requiring HFS have usually simpler 
structures (excepting Cu), and it is safe to conclude that the use 
of HFS has not introduced any important error. 

The atmospheric parameters of the a Cen stars were de- 
termined by simultaneously realizing the excitation & ioniza- 
tion equilibria of Fe I and Fe II. For each star, T e g was ob- 
tained by forcing the Fe I line abundances to be independent 
of their excitation potential. Surface gravity was determined 
by forcing the lines of Fe I and Fe II to yield the same abun- 
dance. Microturbulence velocities £ t were determined forc- 
ing the lines of Fel to be independent of their W\s. The Fe 
abundance [Fe/H] (we use throughout the notation [A/B] = log 
N(A)/N(B) st ar - log N(A)/N(B) Sun , where N denotes the num- 
ber abundance) is automatically obtained as a byproduct of this 
method. The solution thus obtained is unique for a given set of 
gf values, Was and model atmospheres, being independent of 
the starting point and the iteration path. The spectrum synthe- 
sis code is originally due to Dr. Monique Spite (Observatoire de 
Paris-Meudon, Spite 1967), having been continuously up-dated 
in the last 40 years. 

Formal errors are estimated as follows: for T e g, the ler un- 
certainty of the slope of the linear regression in the [Fe/H] vs. 



6 Porto de Mello, Lyra & Keller: Abundance analysis of a Cen A and B 

Table 2. Spectral lines of the elements used in this study. A and B refer to a Cen A and a Cen B. W\S are the raw measurements, 
prior to the correction to the Voigt W\ scale (see text). 



A 


X 


logs/ 




W x (mA) 




A 


X 


logs/ 




W x (mA) 




A 


X 


logs/ 




W x (mA) 




(A) 


(eV) 




Moon 


A 


B 


(A) 


(eV) 




Moon 


A 


B 


(A) 


(eV) 




Moon 


A 


B 






Na 1 








5381.020 


1.57 


-1.855 


65.9 


79.0 


80.9 


5701.557 


2.56 


-2.116 


89.2 


99.5 


127.9 


6154.230 


2.10 


-1.532 


411.6 


62.5 


92.3 


5418.756 


1.58 


-2.116 


52.8 


66.9 


55.8 


5705.473 


4.30 


-1.427 


42.1 


53.0 


66.8 


6160.753 


2.10 


-1.224 


60.5 


86.3 


120.4 






V 1 








5731.761 
5784.666 


4.26 
3.40 


-1.1 15 
-2.487 


60.9 
33.0 


71.1 

42.8 


85.9 
61.(1 






Mg 1 








5657.436 


1.06 


-0.889 


8.5 


10.4 


-10.4 


5811.916 


4.14 


-2.383 


1 1.3 


17.6 


26.3 


5711.095 


4.34 


-1.658 


1 12.2 


124.3 




5668.362 


1.08 


-0.940 


7.3 


1 1.7 


38.2 


5814.805 


4.28 


-1.851 


23.7 


33.6 


44.7 


5785.285 


5.11 


-1.826 


59.(1 


68.3 


88.7 


5670.851 


1.08 


-0.396 


21.6 


31.4 


82.1 


5835.098 


4.26 


-2.085 


16.2 


22.8 


35.(1 














5727.661 


1 .05 


-0.835 


9.8 


14.6 


56.8 


5849.681 


3.69 


-2.963 


8.3 


13.5 


22.9 






Si 1 








6135.370 


1 .05 


-0.674 


14.1 


20.6 


59.6 


5852.222 


4.55 


-1.180 


43.2 


54.3 


71.3 


5517.533 


5.08 


-2.454 


14.5 


25.5 


21.8 


6150.154 


(1.30 


-1.478 


12.9 


20.4 


65.5 


5855.086 


4.61 


-1.521 


24.5 


34.8 


42.8 


5665.563 


4.92 


-1.957 


43.0 


57.0 


62.3 


6274.658 


0.27 


-1.570 


11.1 


15.0 


60.6 


5856.096 


4.29 


-1.553 


36.6 


48.1 


61.2 


5684.484 


4.95 


-1.581 


65.1 


79.4 


78.6 


6285.165 


0.28 


-1.543 


1 1.5 


25.1 


66.3 


5859.596 


4.55 


-0.579 


77.4 


88.4 


105.6 


5690.433 


4.93 


-1.627 


63.2 


67.4 


64.6 














6098.250 


-1.56 


-1.760 


17.7 


25.7 


35.4 


5701.108 


4.93 


-1.967 


41.9 


57.3 


58.2 






Cr I 








6120.249 


0.92 


-5.733 


7.5 


1 1.8 


29.1 


5708.405 


4.95 


-1.326 


78.7 


93.1 


95.9 


5214.144 


3.37 


-0.739 


18.2 


27.2 


39.9 


6137.002 


2.211 


-2.830 


72.9 


83.4 


104.9 


5793.080 


4.93 


-1.896 


-16.1 


62.4 


60.3 


5238.964 


2.71 


-1.312 


19.9 


33.8 


55.2 


6151.616 


2.18 


-3.308 


51.1 


60.7 


811.5 


6125.021 


5.61 


-1.496 


34.7 


51.3 


51.4 


5272.007 


3.45 


-0.311 


32.9 


42.9 


88.0 


6173.340 


2.22 


-2.871 


70.2 


83.7 


100.9 


6142.494 


5.62 


-1.422 


38.5 


52.5 


49.8 


5287.183 


3.44 


-0.822 


13.8 


1 6.6 


34.9 


6219.287 


2.20 


-2.412 


93.6 


102.9 


134.6 


6145.020 


5.61 


-1.397 


-10.5 


55.1 


55.3 


5296.691 


0.98 


-1.343 


96.5 


107.0 


153.1 


6226.730 


3.88 


-2.068 


31.2 


40.1 


56.8 


6243.823 


5.61 


-1.220 


51.8 


67.5 


64.1 


5300.75 1 


0.98 


-2.020 


64.7 


76.6 


103.7 


6240.645 


2.22 


-3.295 


50.0 


60.0 


80.4 


6244.476 


5.61 


-1.264 


48.8 


66.5 


66.3 


5304.183 


3.46 


-0.701 


16.8 


25.1 


16.8 


6265.131 


2.18 


-2.537 


88.3 


99.9 


135.0 














5318.810 


3.44 


-0.647 


19.2 


27.9 


53.4 


6271.283 


3.33 


-2.703 


26.5 


36.7 


54.9 






C.i I 








5784.976 


3.32 


-0.360 


34.0 


48.0 


72.9 














5261.708 


2.52 


-0.564 


126.9 


123.0 


182.1 


5787.965 


3.32 


-0.129 


49.4 


59.8 


83.8 






He 11 








5867.572 


2.93 


-1.566 


26.7 


35.8 


58.6 














5234.630 


3.22 


-2.199 


90.4 


110.9 


92.7 


6161.295 


2.52 


-1.131 


135.5 


85.7 


128.7 






Cr 11 








5264.812 


3.33 


-2.930 


53.1 


67.2 


64.8 


6163.754 


2.52 


-1.079 


126.6 


87.7 


93.8 


5305.855 


3.83 


-2.042 


27.1 


37.8 


24.1 


5325.560 


3.22 


-3.082 


51.1 


66.3 


52.7 


6166.440 


2.52 


-1.116 


72.6 


88.9 


114.0 


5313.526 


4.07 


-1.539 


38.7 


49.3 


43.4 


5414.075 


3.22 


-3.485 


33.7 


46.0 


32.1 


6169.044 


2.52 


-0.718 


97.7 


113.4 


158.0 














5425.257 


3.20 


-3.229 


45.5 


58.1 


46.3 


6169.564 


2.52 


-0.448 


118.5 


137.5 


185.8 






Mn 1 








6149.249 


3.89 


-2.761 


37.8 


-18.4 


30.(1 














5394.670 


0.00 


-2.916 


83.4 


103.8 


165.0 


6247.562 


3.89 


-2.325 


57.1 


69.6 


48.3 






Sc 1 








5399.479 


3.85 


-0.045 


42.9 


63.4 


96.6 














5671.826 


1.45 


0.538 


16.1 


24.2 


64.1 


5413.684 


3.86 


-0.343 


28.0 


45.8 


73.1 






Co I 








6239.408 


0.00 


-1.270 


7.8 


12.3 


12.9 


5420.350 


2.14 


-0.720 


88.6 


1 16.4 


177.5 


5212.691 


3.51 


-0.180 


20.0 


34.8 


50.6 














5432.548 


0.00 


-3.540 


54.9 


73.3 


141.0 


5301.047 


1.71 


-1.864 


22.5 


33.8 


57.9 






Sc 11 








5537.765 


2.19 


-1.748 


37.3 


57.4 


113.3 


5342.708 


4.02 


0.661 


35.1 


48.5 


80.5 


5318.346 


1 .36 


-1.712 


16.0 


23.3 


28.7 














5359.192 


4.15 


0.147 


1 1.9 


20. 1 


37.3 


5526.815 


1.77 


0.099 


80.2 


97.8 


86.5 






Fe I 








5381.772 


4.24 


0.000 


8.2 


15.6 


211.0 


5657.874 


1.51 


-0.353 


71.2 


85.9 


76.5 


5054.647 


3.64 


-1.806 


53.0 


67.6 


97.6 


5454.572 


4.07 


0.319 


18.9 


26.8 


42.8 


5684.189 


1.51 


-0.984 


-10.4 


55.6 


48.9 


5067.162 


4.22 


-0.709 


85.4 


90.3 


123. 1 














6245.660 


1.51 


-1.063 


37.6 


5 I.I 


51.0 


5109.649 


4.30 


-0.609 


87.9 


100.1 


143.0 






M 1 




















5127.359 


0.93 


-3.186 


109.6 


117.7 


164.8 


5094.406 


3.83 


-1.088 


32.4 


-15.0 


62.2 






Ti 1 








5151.971 


1.01 


-3.128 


108.9 


120.8 


177.5 


5220.300 


3.74 


-1.263 


28.5 


411.6 


52.1 


5071.472 


1.46 


-0.683 


36. 1 


45.7 


99.0 


5213.818 


3.94 


-2.752 


7.5 


13.0 


20.2 


5435.866 


1.99 


-2.340 


57.4 


73.8 


88.9 


5113.448 


1.44 


-0.815 


30.9 


36.4 


98.3 


5223.188 


3.63 


-2.244 


32.3 


41.6 


56.5 


5452.860 


3.84 


-1.420 


19.0 


29.1 


37.3 


5145.464 


1.46 


-0.615 


39.6 


50.8 


90.6 


5225.525 


0.1 1 


-4.577 


83.2 


92.7 


128.8 


5846.986 


1.68 


-3.380 


24.2 


35.8 


56.9 


5147.479 


0.00 


-1.973 


43.7 


66.8 


102.7 


5242.491 


3.63 


-1.083 


92.3 


107.5 


132.2 


6176.807 


4.09 


-0.315 


61.2 


86.8 


9(1. 1 


5152.185 


0.02 


-2.130 


34.9 


43.5 


79.1 


5243.773 


4.26 


-0.947 


69.6 


84.3 


91.5 


6177.236 


1.83 


-3.476 


16.4 


29.4 


41.2 


5211.206 


0.84 


-2.063 


9.3 


13.7 




5250.216 


0.12 


-4.668 


78.1 


96. 1 


136.0 














5219.700 


0.02 


-2.264 


28.9 


39.0 


82.9 


5321.109 


4.43 


-1.191 


48.0 


59.5 


76.8 






Cu I 








5295.780 


1.07 


-1.633 


14.2 


18.6 


49.8 


5332.908 


1 .56 


-2.751 


102.8 


118.6 


139.8 


5218.209 


3.82 


0.293 


55.9 


71.8 


82.5 


5426.236 


0.02 


-2.903 


9.0 


12.1 


48.2 


5379.574 


3.69 


-1.542 


64.8 


76.5 


94.4 


5220.086 


3.82 


-0.630 


15.5 


25.0 


32.5 


5679.937 


2.47 


-0.535 


8.6 


10.1 


28.8 


5389.486 


4.41 


-0.533 


87.3 


102.2 


123.0 














5866.452 


1.07 


-0.842 


49.6 


64.8 


107.5 


5395.222 


4.44 


-1.653 


25.2 


33.7 


52.1 






Y □ 








6098.694 


3.06 


-0.095 


6.9 


10.1 


27.6 


5432.946 


4.4-1 


-0.682 


76.4 


90.2 


106.7 


5087.426 


1.08 


-0.329 


49.3 


53.7 


54.2 


6126.224 


1.07 


-1.358 


25.3 


31.9 


68.3 


5491.845 


4.19 


-2.209 


14.3 


25.4 


32.9 


5289.820 


1.03 


-1.847 


4.5 


5.6 


7.7 
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Fig. 3. Left. Fe I line abundances as a function of excitation potential x m a Cen A. The dashed line is the regression of the two 
quantities. By forcing the angular coefficient to be null, we retrieve the excitation T c g of the star. 

Right. Fe I line abundances as a function of equivalent width W\ in a Cen A. A null angular coefficient provides the microturbu- 
lence velocity The ionization equilibrium between Fe I and Fe II constrains the surface gravity, and is realized simultaneously 
with the T c ft and £t determinations. 



X diagram yields the T e g variation which could still be ac- 
cepted at the ler level. For the microturbulence velocity, the 
same procedure provides the ler microturbulence uncertainty 
in the [Fe/H] vs. W\ diagram. For the metallicity [Fe/H], we 
adopt the standard deviation of the distribution of abundances 
derived from the Fel lines. The error in \ogg is estimated by 
evaluating the variation in this parameter which produces a 
disagreement of ler between the abundances of Fe I and Fe II, 
where we regarded the abundance offset as ltr when its value 
was equal to the largest dispersion of the Fe abundances (usu- 
ally that of Fe II). The results of this procedure are shown in 
Fig. [3] where we plot the iron abundances of a Cen A derived 
from lines of both Fe I and Fe II against the excitation potential 
and W\s. The baseline of the Fe I lines is seen to be large both 
in W\ and \. 

Additional effective temperatures were determined by fit- 
ting the observed wings of Ha, using the automated proce- 
dure described in detail in Lyra & Porto de Mello (2005). We 
have employed for the a Cen stars new spectroscopic data from 
the Observatorio do Pico dos Dias, with the same resolution 
but greater signal-to-noise ratio than used by Lyra & Porto de 
Mello (2005). This procedure is shown in Fig. [4] and Fig. [5] 
We found T cB = 5793 ± 25 K (a Cen A) and T cS = 5155 ± 4 
K (a Cen B). The moon spectrum is very well fitted by the pa- 
rameters adopted for the NMARCS solar atmosphere model. 
The quoted standard errors refer exclusively to the dispersion 
of temperature values attributed to the fitted profile data points. 
This makes the uncertainty of the a Cen B T c g artificially very 
low, due to the high number of rejected points. An analysis of 
errors incurred by the atmospheric parameters assumed in the 
fitting procedure, plus the photon statistics (not including pos- 
sible systematic effects produced by the modelling, see Lyra & 
Porto de Mello (2005) for a full discussion), points to an aver- 
age error of —50 K in the T e gS determined from the Ha line. 



For the very low-noise spectra of these two stars, the expected 
errors would be slightly less, but the greater difficulty in find- 
ing line-free sections in the Ha profile of the severely blended 
spectrum of a Cen B offsets this advantage. For the latter, thus, 
the probable uncertainty should be closer to ~ 100 K. The nor- 
malization procedure is a a relevant source of error for T c g& de- 
rived from Ha, as discussed by Lyra & Porto de Mello (2005): 
they found that a 0.2% error in the continuum level translates 
to -25 Kin T cff . 

As an external check on our normalization procedure, we 
have compared Ha spectra of the Sun (moon) from Lyra & 
Porto de Mello (2005), from eight independent observing runs, 
also reduced independently, with the data of ABLC04, who 
performed a careful two-dimensional continuum normalization 
of the echelle spectra. Due to normalization problems in the 
FEROS/ESO (La Silla) spectra, these authors employed Ha 
spectra as a T e s criterion only for the northern stars of their 
sample, which could be observed with the McDonald 2dcoude 
spectrograph. So a direct comparison between the two sets of 
Ha spectra is not possible for the a Cen stars. We found, for the 
solar spectra, an average difference of only (—0.29 ± 0.39) %, 
as measured in those regions relevant to the T e s determination. 
This assures us of the absence of important systematic errors 
in this respect. A comparison of the Ha normalized spectra is 
shown in Fig. [6] 

Another check on the T e g values of the two stars may 
be obtained from the IRFM scale of Ramirez & Melendez 
(2005a). Ramirez & Melendez (2005b) compare, for the two 
a Cen stars, direct T e fjs, obtained from measured bolometric 
fluxes and angular diameters, with those determined from the 
IRFM method, as well those obtained from the application of 
their own T c ff(IRFM) scale to color indices and an adopted 
metallicity of [Fe/H] = +0.20, for both stars. Respectively, they 
find, for a Cen A, T cff (direct) = 5771, T off (IRFM) = 5759 K 
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and T ff(calibration) = 5736 K; the corresponding values for 
aCenB are, respectively, T eS (direct) = 5178 K, T off (IRFM) = 
5221 K and T ff(calibration) = 5103 K. They adopt as weighted 
averages T eff (a Cen A) = 5744 ± 72 K and T cff (aCenB) = 
5136±68K. For a Cen A, a good agreement is found between 
this T ff and our Ha one. Formally, there is also a reasonable 
agreement between this T c g and our spectroscopic one. For 
a Cen B, however, the spectroscopic T c g is significantly higher 
than those derived from Ha and the IRFM method. 

Yildiz (2007) drew attention to often neglected BVRI 
(Cousins system) measurements of the system's components 
by Bessell (1990). It has long been considered risky to use 
color indices for the T c g determination of very bright stars, in 
which a variety of systematic effects are expected as compared 
to standard stars of photometric systems, among which non- 
linearity, detector dead-time, and, in the case of a Cen, possible 
contamination by the companion (see Chmielewski et al. 1992 
for a full discussion). Introducing Bessell's color indices into 
the Ramirez & Melendez (2005b) calibrations, along with our 
metallicities (Tabled, we derive T cff (a Cen A) = 5794 ± 34 K 
and T ff(aCenB) = 5182 ± 19 K, as a weighted average of the 
(B — V), (V — R) and (V — I) color indices, the latter two in the 
Cousins system. These new photometric T e g determinations, 
directly from the calibrations, agree well, for a Cen A, with 
both the Ha and the spectroscopic one. On the other hand, for 
a Cen B this determination lessens slightly, but does not elim- 
inate, the disagreement between the spectroscopic T g and the 
other two. We must therefore state clearly that there is an offset 
between the spectroscopic T g of a Cen B and the other two 
T e ff determinations. These last figures, in a classical spectro- 
scopic analysis of solar-type stars, would be regarded as "pho- 
tometric" ToffS, to be compared to those obtained from other 
methods. Our results are all displayed in Table [3] where we 
also list the direct surface gravities resulting from the observed 
radii and dynamical masses. 

Valenti & Fischer (2005) have derived the atmospheric pa- 
rameters of a Cen A and B by means of a different technique. 
They have fitted directly large sections of the observed spectra 
to synthetic ones, obtaining the atmospheric parameters (this 
technique also relies on the excitation & ionization equilibria 
of atomic species). Their analysis is differential with respect to 
the Sun, for which they adopted T cS = 5770, [Fe/H] = +0.00 
and log g = 4.44. They quote uncertainties of 44 K, 0.03 dex 
and 0.06 dex, respectively, for T e g, log g and [Fe/H]. Their 
results are T cff = 5802 K, log g = 4.33 and [Fe/H] = +0.23, for 
a Cen A, and T off = 5 178 K, log g = 4.56 and [Fe/H] = +0.21, for 
a Cen B. These figures are in good agreement, even within their 
very small claimed uncertainties, with our spectroscopic pa- 
rameters, again, with the exception of the spectroscopic T c g of 
aCenB. Particularly, their log g values are in excellent agree- 
ment with the direct log g values of TableO 

Frutiger et al. (2005) have also analysed the spectra of 
a Cen A and B by a fundamentally different and promising 
technique. They have inverted high resolution (R ~ 10 5 ), mod- 
erately high S/N (~250) spectra of the stars by means of 
a multi-component model photosphere. The components take 
into account rotational broadening, center-to-limb variations 
and vertical and horizontal flows of surface elements, such as 



granules and inter-granular areas. In this approach, the full line 
profile is used to constrain the temperature stratification of the 
atmosphere, as well as the velocity fields (Allende-Prieto et 
al. 1998). The technique is rather model-dependent, however, 
and should be compared with classical spectroscopic analyses 
with caution. For the 3-component models, the ones they fa- 
vor, T ffS and [Fe/H] substantially lower than those found by 
us are obtained (their table 4). For a Cen A, they favor T c g= 
5705 K, log g = 4.28 ± 0.03 and [Fe/H] = +0.08 ± 0.02 (we 
have converted their abundances from absolute to relative val- 
ues with respect to the Sun, which they also analysed with the 
same techniques. Their analysis, in this sense, may also be con- 
sidered as differential). For aCenB, the results are T c g= 5310 
K, log g = 4.74 ± 0.02 and [Fe/H] = +0.05 ± 0.01. 

It is not straightforward to determine the uncertainties of 
their T c g values, since they only quote the uncertainty of 
this parameter as derived by a weighted average of the ctT 4 
fluxes of each atmospheric component, weighted by the filling- 
factors, by means of the Eddington-Barbier relation applied 
to a grey atmosphere. The T e fjS we quoted are the ones they 
computed from the flux spectra obtained from their temper- 
ature stratifications fed into an ATLAS9 model, for which 
they provided no formal uncertainty. The uncertainties of their 
T ffS within the Eddington-Barbier approximation, are, respec- 
tively, 130 K and 67 K for a Cen A and B. Within these error 
bars, then, their T c gs resulting purely from the inversion proce- 
dure (letting all parameters free) may be regarded as compat- 
ible with ours. Their metallicities, however, are significantly 
lower. They also found a very high log g for a Cen B, prompt- 
ing them to attempt a larger number of inversions this for com- 
ponent, first fixing log g = 4.48 dex, which led to a reduced 
T c ft= 5154 K with no appreciable change in metallicity, and 
then fixing the rotational velocity, which produced T c g= 5260 
K, again with no significant impact on metallicity, and a new 
log g = 4.68 ± 0.05. 

The T e gS favored by the inversion method of Frutiger et 
al. (2005) add to a complex situation. They seem to be in good 
agreement with our spectroscopic T c g for a Cen B when all pa- 
rameters are independently derived from the inversion method, 
but the agreement switches to one with the Ha and photometric 
T c ff s, when the surface gravity is fixed. All [Fe/H] values they 
obtain are lower than ours. For a Cen A, their surface gravity 
is in good agreement with ours and the direct one. It is dif- 
ficult, however, to reconcile their surface gravity for a Cen B 
with our spectroscopic one, and the direct one given in Table[3] 
It is not clear whether their results can be directly compared to 
ours, given the difference in approach. These authors discuss 
the possibility of improving their technique by higher resolu- 
tion spectra: they remark on the difficulty of disentangling ef- 
fects of rotation, macro-turbulence, granulation and instrumen- 
tal profile. This sophisticated approach, probably, can lead to 
substantially improved constraints to the atmospheric parame- 
ters of solar-type stars, leading to increased physical insight on 
the shortcomings of 1-D, static atmospheric models. 
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Fig. 4. Effective temperature determination by fitting theoretical profiles to the wings of Ha for a Cen A. Upper panel: crosses 
refer to pixels eliminated by the statistical test and 2a criterion (see Lyra & Porto de Mello 2005 for details). The T e g derived by 
the accepted pixels (open circles) is 5793 K. The corresponding best determined line profile is overplotted (solid thick line) on 
the observed spectrum. The dashed lines are theoretical profiles, spaced by 50 K and centered in 5847 K, the spectroscopic T e g . 
Lower panel: the spectrum at a larger scale. 



3.1. Systematic offsets between spectroscopic, Balmer 
line and photometric T c gscales 

The study of the atmospheric parameters collected at Table Q] 
reveals an interesting pattern, if one again considers only the 
analyses published since the 90s, invariably based on high S/N 
spectra acquired with solid-state detectors. For a Cen A, the 
works employing the Fe i/Fe II criterion, alone or combined 
with another method, obtained T e gs generally higher than those 
derived exclusively from photometry (excepting the result of 
FM90, which is one of the lowest and applied the excitation 
& ionization approach). Thus, the present work, Santos et al. 
(2004) and Neuforge-Verheecke & Magain (1997) found the 
highest T ffS, while ABLC04, relying only in photometry, de- 
rived the lowest T c g . Doyle et al. (2005) also used the stellar lu- 
minosity and radius to derive T c g , and found a value lower than 
the spectroscopic ones, but for the FM90 analysis. ABLC04 
reported good agreement of their photometric T c g s with those 
derived from the H/3 profiles. Interestingly, they also found an 
offset of ^120 K between their T c gs and the ones of Kovtyukh 
et al. (2003), the latter being the higher. Kovtyukh et al. (2003) 
derived T c g s by the line-depth ration method (Gray & Johanson 



1991), and their T c ffS were actually calibrated by spectroscopic 
ones. The pattern just discussed suggests that the spectroscopic 
T e g scale is indeed hotter than the photometric and Balmer line 
one, and that the latter two are in generally good agreement. 
The situation for aCenB is unfortunately much less clear: 
there are less analyses, and the authors employed a more re- 
stricted set of criteria. Nevertheless, the lowest T ff is again 
due to photometric methods (ABLC04), the highest one in this 
case corresponding to the Ha-derived T Q g of Chmielewski et 
al. (1992). 

A disagreement between photometric and spectroscopic 
T e g scales has been recently pointed out by a number of au- 
thors. Ramirez et al. (2007) discussed how the Fe i/Fe II ion- 
ization equilibrium is not realized in cool stars when an IRFM 
T e g scale is applied, in the derivation of oxygen abundances 
from the A7774 triplet lines. Their sample is large, and they 
convincingly show (their fig. 5) that the T c g offset between 
the IRFM and the spectroscopic scales is significant for T c g ~ 
5000 K (a reasonable agreement was found for T c g ~ 6000 K). 
They note that Ramirez & Melendez (2004) found the same 
offset, as did Santos et al. (2004), Yong et al. (2004) and Heiter 
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Fig.5. The same as in Fig|4]for aCenB. The T c ff derived by the accepted pixels is 5155 K. The corresponding line profile is 
overplotted (solid thick line) on the observed spectrum. The dashed lines are theoretical profiles, spaced by 50 K and centered in 
53 16 K, the spectroscopic T e ff ■ Lower panel: the spectrum at a larger scale. 



& Luck (2003). Yong et al. (2004) suggested that non-LTE ef- 
fects, shortcomings of the model atmosphere representation for 
cool stars, or as yet unidentified effects, might be responsible 
for the discrepancy. Ramirez et al. (2007), however, note that 
Santos et al. (2005) reported good agreement between IRFM 
and spectroscopic T c ff scales. Adding to this complex picture, 
Casagrande et al. (2006), in the derivation of their own IRFM 
T e ff scale, also found good agreement between spectroscopic 
and IRFM T c ff s. They argue that the disagreement reported by 
other authors might be due, at least in part, to uncertainties 
in the different absolute flux calibrations adopted, and suggest 
that additional direct angular diameter measurements for a well 
chosen sample of G and K dwarfs might go a long way towards 
clarifying the disagreement of the T ff scales. 

The results of Casagrande et al. (2006) were essentially 
backed by the T c ff scale of Masana et al. (2006). The latter 
employed a variation of the IRFM method, in which the stellar 
energy distributions were fitted, from the optical to the IR, to 
synthetic photometry computed from stellar atmosphere mod- 
els. They found their results only slightly offset from the IRFM 
results. The offset, at ~30 K, was deemed to be small, and these 
authors optimistically assert that their T e g scale agrees with 
the spectroscopic one for FGK dwarfs and subgiants, such that 



T ffS for this class of stars may be regarded as accurate within 
— 1% or better. 

The consistency of the different T c ff scales is sought as an 
important confirmation that 1-D, plane-parallel, static and LTE 
model atmospheres adequately represent cool stars, if not in the 
absolute, at least in the relative sense, provided that the Sun can 
be accurately placed in the stellar context. The solar placement 
in the Fe i/Fe II excitation & ionization equilibria and Balmer 
line T c ff scale is obtained by the observation of solar flux spec- 
tra. An accurate photometric placement of the Sun in the cor- 
responding T c ff scale, however, is a difficult task still beset 
with large uncertainties (see, e.g., Holmberg et al. 2006 for an 
up-to-date discussion). These three T e ff scales actually gauge 
rather dissimilar physical quantities. The excitation & ioniza- 
tion T G ff is obtained by matching models to observed spectral 
line intensities. The Balmer line T c ff measures the temperature 
stratification of the atmosphere, which is mapped onto the line 
wings by the depth-dependence of the source-function. The 
photometric T c ff must reproduce the stellar flux distribution in 
a large wavelength regime, and is the one most directly tied 
to the fundamental definition of effective temperature (Bohm- 
Vitense 1981). As long as the consistency between these scales 
is realized no better than within ^ 150 K, the T G ff of cool dwarfs 
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Fig. 6. A comparison of the normalization of our Ha spectrum of the Sun (moon) with that of ABLC04. The large differences 
in the telluric lines are apparent. The left and right panels depict the regions outside the Ha profile which were used for the 
continuum normalization. In the center panel we depict, as the shaded gray areas, the wavelengths of the wing profile used for 
the T e ff determination. Note that the abscissa is truncated to show only the most relevant portions of the spectra, the blue and red 
limits, which set the normalization scale, and those used directly for the fit of the line wings. In this example, the mean difference 
between the two spectra (computed only for the line wing regions actually fitted) is -0.8 ± 0.7 %, and it is the worst case of our 
comparison of eight spectra (see text for details). 



and subgiants will remain uncertain by this amount at the very 
least. 

Non-LTE and other possibly more complex effects have 
repeatedly been blamed for offsets between spectroscopically 
and photometrically derived atmospheric parameters in cool 
stars. Schuler et al (2006a), in their analysis of Hyades dwarfs, 
reported a systematic offset of the oxygen abundances derived 
from the A7774 triplet lines, for T e s < 5450 K, in the oppo- 
site sense of the NLTE expectations. They tentatively suggest 
that chromospheric activity might be at least partially respon- 
sible for the offset, an explanation also concurred by Morell 
& Micela (2004), though the latter propose that model atmo- 
sphere pitfalls might be also present. Schuler et al. (2006b) 
reinforce this interpretation in an analysis of the A6300 [OI] 
line in the very active Hyades stars, reporting offsets between 
Fe i/Fe II abundances which increase as T c ff < 5000 K. It 
should be emphasized that high chromospheric activity is un- 
likely to be a source of the T G ff discrepancy of a Cen B, since 
both components of the system are inactive stars, which prob- 
ably implies that the problem is more complex. We also draw 
attention to the result of Shchukina & Trujillo-Bueno (2001), 
who found an offset of the Fe i/Fe II abundances of the Sun. 
They interpret this offset as well explained by NLTE effects 
amounting to 0.07 dex for Fe I, the best fit Fe abundance for a 
LTE analysis being the lower by this amount. They assert that a 
full 3-D, NLTE model atmosphere formulation is able to bring 
the solar photospheric Fe abundance in line with the meteoritic 
one, at log N(Fe) = 7.50 (in the usual scale where log N(H) 
= 12.00). The main cause of the offset is the overionization of 
Fe I, and the larger errors are seen in the Was of low excitation 
lines, which are weaker in the NLTE case. 



It is interesting to note that this effect, in a classical LTE 
model atmosphere analysis, would result in the overabundance 
of the high excitation Fe I lines, an effect naturally interpreted 
as too low a T ff being attributed to the model. This is ex- 
actly the condition necessary for a 1-D, LTE analysis to lead 
to the high spectroscopic T c ff — that we obtained. Forcing the 
Fe i/Fe II abundances into agreement in a LTE analysis would 
indeed call for a higher T e g , to a first approximation, by ~ 100 
K, a value similar to the difference between our spectroscopic 
and Ha/photometric T Q ss of a Cen B. We suggest that the simi- 
larity in atmospheric parameters between a Cen A and the Sun, 
in the context of a differential analysis, led to a good agreement 
between the three T c g criteria for the former. For aCenB, a 
much cooler object, an imperfect cancellation of the presence 
of NLTE effects is probably the reason why the three different 
T ff criteria do not agree. 

Even if the presence of NLTE effects and other problems 
can be established, one must keep in mind that other uncertain- 
ties are present in the photometric and Ha T ff scales, as dis- 
cussed above. It would be very valuable to extend the novel ap- 
proaches of Valenti & Fischer (2005) and Frutiger et al. (2005) 
to additional stars for which stringent observational constraints 
are available, to quantify how discrepancies of the type dis- 
cussed here do occur for objects with atmospheric parame- 
ters increasingly different from the Sun's. At present, we can 
state that it is likely that our spectroscopic T e fj for aCenB 
is more uncertain that the other two determinations, being too 
high. This conclusion is, however, drawn in the context of a 
confusing picture. Additional work is clearly necessary before 
a definitive judgment can be passed on the consistency be- 
tween the excitation & ionization, B aimer line and photometric 
T e ff scales can be reached. 
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Table 4. The abundance pattern of the a Cen stars. The first 
two columns give the [X/Fe] ratios with the corresponding ob- 
served dispersions of the line abundances, for those elements 
with three or more available lines, for a Cen A. The next two 
columns give the corresponding data for a Cen B. the fifth col- 
umn provides the total uncertainty of the [X/Fe] ratios corre- 
sponding to errors in T e g, log g, [Fe/H], £ and W\, respec- 
tively, of 30 K, 0.12 dex, 0.04 dex, 0.03 kms" 1 and 2.9 mA, 
composed in quadrature (the latter enters twice, see section 3). 



element 


a Cen A 


aCenB 


A 


[X/Fe] 


a 


[X/Fe] 


a 


Na 


+0.14 




+0.25 




0.03 


Mg 


0.00 




+0.04 




0.03 


Si 


0.00 


0.06 


0.00 


0.05 


0.05 


Ca 


+0.03 


0.03 


+0.09 


0.08 


0.03 


Sc 


0.00 


0.02 


+0.01 


0.07 


0.06 


Ti 


-0.01 


0.04 


+0.19 


0.14 


0.06 


V 


+0.02 


0.05 


+0.21 


0.09 


0.07 


Cr 


0.00 


0.06 


+0.06 


0.09 


0.06 


Mn 


+0.07 


0.01 


+0.19 


0.06 


0.04 


Co 


0.00 


0.06 


+0.06 


0.07 


0.07 


Ni 


+0.10 


0.07 


+0.11 


0.05 


0.04 


Cu 


+0.10 




+0.10 




0.04 


Y 


-0.09 


0.05 


-0.09 


0.08 


0.07 


Ba 


-0.13 


0.05 


-0.16 


0.08 


0.06 



Notwithstanding the discrepancies, if we follow an usual 
practice in spectroscopic abundance analysis which obtain at- 
mospheric parameters with more than one criterion, mean 
T c ffS for the a Cen stars can be calculated from the values of 
Table [3] averaged by their inverse variances. The results are 
T cS (a Cen A) = 5824 ± 27 K and T cS (a Cen B) = 5223 ± 62 
K, where the quoted uncertainties are the standard deviations 
of the average, and do not reflect external and systematic er- 
rors. Good agreement between our ionization surface gravities, 
and those directly determined from measured masses and radii, 
within the errors, for both stars, is realized. 

4. Abundance pattern 

The abundances of the other elements were obtained with the 
adopted atmospheric model of each star, corresponding to the 
spectroscopic T e jps, and the log g, [Fe/H] and £ values as given 
in Table [3] Average abundances were calculated by the straight 
mean of the individual line abundances. For Sc, Ti and Cr, good 
agreement, within the errors, was obtained for the abundances 
of the neutral and singly ionized species, and thus these species 
confirm the Fe i/Fe II ionization equilibria. The results are plot- 
ted in Fig. [7] as [X/Fe] relative to the Sun. In this figure the 
error bars are merely the internal dispersion of the individual 
line abundances. For Na I, Mg I, Cu I and Ba II, the dispersions 
refer to the difference between the abundances of the two avail- 
able lines for each element. There is a very good consistency 
between the abundance patterns of the two stars, but for the 
the Ti I and V I abundances. The larger error bars for a Cen B, 



however, lead us to consider the two abundance patterns as con- 
sistent with each other. 

In Fig. [7] the vertical dark grey bars besides the data points 
of the abundance pattern of a Cen A are the composed r.m.s. 
uncertainties, for each element, calculated by varying the spec- 
troscopic atmospheric parameters of a Cen A by the corre- 
sponding uncertainties of Table[3] To this calculation we added 
the abundance variations caused by summing to all W\& the 2.9 
mA uncertainty of the correction of Fig|2] This W\ uncertainty 
enters twice: once for the uncertainty in the corrected moon 
W\S, reflecting onto the solar log gfs, and another one due to 
stellar W\ themselves. In Fig [7j it is apparent that the abun- 
dance variations due to the uncertainties in the atmospheric pa- 
rameters and W\s are comparable to the observed dispersions 
of the line abundances for a Cen A. For a Cen B, the line abun- 
dance dispersions are generally larger, probably due to its more 
uncertain W\s, but also, as discussed above, possibly owing to 
its larger T e g difference from the standard object (the Sun) and 
its more uncertain atmospheric parameters. 

Our abundance pattern for a Cen A is clearly the most re- 
liable of the pair, and is directly compared to those of other 
authors in Fig [8] Only abundances represented by more than 
one spectral line are shown. The observed dispersion is com- 
parable to the uncertainties normally quoted in spectroscopic 
analyses. Only for the light elements between Mg and Ti is a 
larger disagreement observed, in this case due to the analysis 
of ABLC04, in which abundances are higher than in the bulk 
of other data by ^0.2 dex. For the elements heavier than V, es- 
sentially all data agree that V and Cr have normal abundance 
ratios, that Mn, Co, Ni and Cu are enhanced, and all heavy el- 
ements from Y to Eu are deficient in the abundance pattern of 
a Cen A with respect to the Sun, with the sole exception of Ba, 
for which ABLC04 found a normal abundance. The available 
literature data also suggests that the C, N and O abundance ra- 
tios of a Cen A are solar. 

This statistical analysis can be extended if we regard only 
the elements for which at least three independent studies pro- 
vided data. This is shown in Fig [9] for a more select sample 
of elements. We may conclude, with somewhat greater robust- 
ness, considering the number of abundance results, that Na, 
Mg, Si, Mn, Co and Ni are over-abundant; that Ca, Sc, Ti, V 
and Cr have solar abundance ratios; and that Y and Ba are over- 
deficient in the abundance pattern of a Cen with respect to the 
Sun. 

The high metallicity of the a Cen system, and its space ve- 
locity components (U,V,W)(km.s _1 ) = (-24, +10, +8) (Porto 
de Mello et al. 2006, all with respect to the Sun) place it un- 
ambiguously as a thin disk star. We next analyze its abundance 
ratios, for the elements with more reliable data, as compared to 
recent literature results for metal-rich stars. Bensby et al. (2003, 
hereafter BFL) analysed 66 stars belonging to the thin and thick 
disks of the Milky Way, deriving abundances of Na, Mg, Al, 
Si, Ca, Ti, Cr, Fe, Ni and Zn. Bodaghee et al. (2003, here- 
after BSIM) analysed a sample of 119 stars, of which 77 are 
known to harbor planetary companions, deriving abundances 
of Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co and Ni. The data of the 
latter study comes from the Geneva observatory planet-search 
campaign (e.g., Santos et al. 2005). Since they concluded that 
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Table 3. The atmospheric parameters resulting from our spectroscopic analysis: from the excitation & ionization equilibria of 
Fe I and Fe II (excitation T e fj, ionization log g, [Fe/H] and £ t ), from the wings of Ha and the photometric calibrations of Ramirez 
& Melendez (2005b)(photometric). The mean T e g is weighted by the inverse variances. The "direct" log g values are derived 
from the masses of Pourbaix et al. (2002) and the radii of Kervella et al. (2003) (see text). 





Teff(K) 

excitation 


T eB (K) 
Ha 


Teff(K) 

photometric 


Teff(K) 

weighted mean 


log 9 
ionization 


logs 
direct 


[Fe/H] 


CiCkms" 1 ) 


a Cen A 
aCenB 


5847±2? 
5316±28 


5793±50 
5155±ioo 


5794±34 
5182±is 


5824±26 
5223±62 


4.34±o.i2 
4.44±o.i5 


4.307±o.oos 
4.538±o.oo8 


+0.24±o.o3 
+0.25±o.o4 


1.46±o.o3 
1.28±o.i2 




Fig. 7. The abundance pattern of a Cen A & B. The uncertainty bars are the dispersions of the abundances given by the element 
lines. The gray bars beside the data points in the a Cen A plot correspond to the total compounded errors arising from the 
atmospheric parameters and the Was (see TableE}. Na is seen to be overabundant, while a solar pattern is seen from Mg to Co, 
but for an excess of Mn. Ni and Cu are also overabundant. Some doubt can be cast about Ti and V, since they seem overabundant 
in a Cen B although solar in a Cen A. The slow neutron capture elements Y and Ba are in clear deficit. The bigger uncertainty 
bars seen in a Cen B are probably a result of a less accurate normalization of the spectra of a cooler star, and may also be due to 
its less accurate atmospheric parameters. 



planet-bearing stars merely represent the high metallicity ex- 
tension of the abundance distribution of nearby stars, their full 
sample can be used to adequately represent the abundance ra- 
tios of metal rich stars, without distinction to the presence or 
absence of low mass companions. These two works have in 
common the important fact that they sample well the metallic- 
ity interval +0.20 < [Fe/H] < +0.40, an essential feature for 
our aim. 

The elements in common between the two abundance sets 
are Si, Ca, Ti, Cr and Ni. It can be concluded (Fig. 13 of BFL, 
Fig. 2 of BSIM) that, in [Fe/H] > +0.20 stars, Ca is under- 
abundant, Ti is normal and Ni is enhanced. For Si, BFL suggest 
abundance ratios higher than solar, while BSIM found a normal 
abundance. For Cr, the data of BSIM suggests abundance ratios 
lower than solar, while BFL found solar ratios. For the elements 
not in common in the two studies, Na, Mg, Al, Sc, V, Mn and 
Co are found to be enhanced in [Fe/H] > +0.20 stars, while Zn 
has normal abundance ratios. 

An interesting feature of the a Cen abundance pattern is the 
under-abundance in the elements heavier than Y, which could 
be reliably established for Y and Ba (Fig[9]). This result is con- 
firmed by Bensby et al. (2005), who found both for Y and Ba 



lower than solar abundance ratios for metal rich thin disk stars. 
Another interesting feature, the excess of Cu (found by us and 
FM90), can be checked with the recent results of Ecuvillon et 
al. (2004), who also obtain for [Fe/H] > +0.20 stars an average 
[Cu/Fe] -0.1 dex. 

Merging our evaluation of the abundance pattern of 
a Cen A, from the available independent analyses, with the pre- 
vious discussion, we conclude that a Cen A is a normal metal- 
rich star in its Na, Mg, Si, Ti, V, Cr, Mn, Co and Ni abun- 
dances. The result for Ca is inconclusive, and only for Sc does 
its abundance diverge from the BFL/BSIM data, in that its nor- 
mal abundance ratio contrasts with the overabundance found 
for [Fe/H] > +0.20 stars by BSIM. It seems reasonably well 
established then that the a Cen system is composed of two nor- 
mal metal rich stars when regarded in the local disk population. 

5. Evolutionary state 

An important outcome of the present analysis is to establish 
if the derived atmospheric parameters, coupled to high-quality 
parallaxes, allow a consistent determination of the masses and 
ages of the a Cen system in a traditional HR diagram anal- 
ysis, matching the stringent constraints posed by the orbital 
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Fig. 8. A comparison of the [X/Fe] abundance pattern of a Cen A derived by different authors. E88 stands for Edvardsson (1998). 
Only elements for which abundances are based on more than one line are included. Apart from the ABLC04 results, which point 
to excesses of nearly all elements with respect to the Sun, most published data imply an abundance pattern in which C, N, O, Ca, 
Sc, V and Cr are normal; Na, Mn, Co, Ni and Cu are enhanced; Mg, Al and Si are possibly enhanced; and Ti and all elements 
heavier than Y are under-abundant. 



solution and seismological data. In Fig. [T0]we plot the posi- 
tion of the a Cen components in the theoretical HR diagram of 
Kim et al. (2002) and Yi et al. (2003), corresponding to its ex- 
act metallicity and a solar abundance pattern. The bolometric 
corrections were taken from Flower (1996), and in calculating 
the luminosities the Hipparcos (ESA 1997) parallaxes and vi- 
sual magnitudes were used. Evolutionary tracks and isochrones 
from different authors (e.g., Girardi et al 2000, Charbonnel et 
al. 1999, Schaller et al. 1992) were also tested, and good agree- 
ment between the different tracks was found, to better than ^50 
K, for the position of both a Cen components. This is an ex- 
pected result, given that the Sun is generally used to calibrate 
these models. The solar mass, radius and age provides a zero 
point to the models and allow a solution as a function of the 
adopted mixing-length of the convection theory (still a free pa- 
rameter) and the initial helium abundance. Thus, differences 
between the models can be substantial in the treatment of stars 
which are very different from the Sun (e.g., Lyra et al. 2006), 
but good agreement for solar-type stars is a natural outcome of 
this procedure. The conclusions drawn below, then, are essen- 
tially model-independent, at least for the T c g and luminosity 
intervals involved here. 

The T c fi values and error bars in the diagram are those of 
the weighted mean of table[3] From the diagram, masses of Ma 
= 1.13 ± 0.01 and M B = 0.89 ± 0.03 can be derived, and agree 
well with the orbital solution of Pourbaix et al. (2002). The age 
of a Cen A can be relatively well constrained to the interval of 
4.5 to 5.3 Gyr {la), again, in good agreement with the seis- 
mological results of Yildiz (2007), Eggenberger et al. (2004), 
Miglio & Montalban (2005) and Thoul et al. (2003), within the 



quoted uncertainties. We conclude that, adopting the average 
T c ffS and [Fe/H] found in this work, the position of a Cen A in 
up-to-date theoretical HR diagrams can be reconciled both with 
seismological and dynamical data. Despite its higher mass and 
its being (probably) older than the Sun, the higher metallicity 
slows the evolution to the point that the star has not yet reached 
the "hook" zone of the HR diagram, thus enabling a unique age 
solution through this type of analysis. 

Concerning a Cen B, it is also apparent in Fig [10] that its 
position cannot be reconciled, within la, with an age near 5 
Gyr, as was possible for a Cen A. However, an upward re- 
vision of only ^60 K would bring its position in agreement 
with a track of 0.93 solar mass (the dynamical value), and an 
age of ~5-8 Gyr. Given the uncertainties discussed in section 
3, along with the probable effect of systematic errors in the 
T e g determination of a Cen B, this is not outside the 2a con- 
fidence interval of the results. We conclude that a Cen A has 
its position in the theoretical HR diagram well matched by up- 
to-date models, within the uncertainties of the determination 
of its atmospheric parameters, and also within the small dif- 
ferences, in this T ff and luminosity regime, between differ- 
ent grids of evolutionary models, a Cen B, in its turn, has been 
shown to suffer from a larger uncertainty in its T c g, precluding 
a more stringent assessment of a match between its evolution- 
ary mass and age and the results from the dynamical solution 
and asteroseismology. Further data are still necessary to allow 
a more definitive conclusion on this issue. If indeed its spectro- 
scopic T c ft is systematically offset, the mean T c g we derived 
would decrease and displace its position on the HR diagram 
to the right, forestalling a match with the age of a Cen A. For 



Porto de Mello, Lyra & Keller: Abundance analysis of a Cen A and B 

a Cen A 



15 




N Na Al Ca Ti Cr Co Cu Zr Ce Elf 



O Mg Si Sc V Mn Ni 



Ba Sm 



Fig. 9. The same as Fig(9j but averaging the [X/Fe] abundances of the elements with published data from at least three authors. 
The boxes comprise 50% of the data points, and are centered at the mean. The horizontal dashes inside the boxes mark the 
median. The whiskers mark the maximum and minimum values. In this more stringent statistical analysis of the available data, 
Na, Mg, Si, Mn, Co and Ni are enhanced; Ca, Sc, Ti, V and Cr are normal; Y and Ba are under-abundant. 



a better understanding of the onset of possible NLTE effects 
in cool stars, and the hindrance thus incurred in the determi- 
nation of their atmospheric parameters, it would be interesting 
to perform further analyses of such objects, for which high- 
quality spectra could be secured, and the relevant observational 
constraints made at least partially available. Interesting bright, 
nearby K-dwarf candidates for such an enterprise are e Eri, 36 
and 70 Oph, o 2 Eri and a Dra. 

6. Conclusions 

We have undertaken a new detailed spectroscopic analysis of 
the two components of the a Centauri binary system, and have 
attempted an appraisal of the many discordant determinations 
of its atmospheric parameters and abundance pattern, and of 
the sources of errors in their determination. We derived purely 
spectroscopic atmospheric parameters, from R = 35 000 and 
S/N > 1 000 spectra, in a strictly differential analysis with the 
Sun as the standard. We obtained T c sa = 5847 ± 27 K and 
T e s b = 53 16 ± 28 K from the spectroscopic analysis, and T e ff a 
= 5824 ± 26 K and T cSB = 5223 ± 62 K from the average of 
spectroscopic, Ha and photometric T e ffs). We derived [Fe/H] 
= +0.24 ± 0.03 dex for the system. The spectroscopic surface 
gravities, log g a = 4.34 ± 0.12 and log g b = 4.44 ± 0.15, are a 
good match to those determined from directly measured masses 
and radii. Good agreement, in both components, is found be- 
tween the photometric T c g and the one resulting from the fit- 
ting of the wings of Ha. For a Cen A, these two T c ff s also agree 
with the spectroscopic one. However, for a Cen B, the spectro- 
scopic T e s was found to be significantly higher, by ^140 K, 
than the other two. A comparison of the published T c ffS for 



the system in the last 20 years roughly support a spectroscopic 
T e g scale hotter than the ones owed to photometric methods or 
the fitting of Balmer lines. 

A comparison with recent results from other techniques 
revealed an unclear picture. Atmospheric parameters for the 
a Cen stars derived by Valenti & Fisher (2005) by fitting di- 
rectly synthetic spectra to a large spectral coverage, agree 
well with our determinations, but for the spectroscopic T c g of 
aCenB. Their surface gravities and metallicities are also in 
line with our figures. Frutiger et al. (2005), inverting high- 
resolution line profiles, found a substantially lower [Fe/H], and 
their model-dependent T c g s agree either with our spectroscopic 
or with the photometric/Ha T c g , depending on assumptions. 
Their log g for a Cen B is also higher than all other recent de- 
terminations. 

We discuss possible origins of the offset between the 
T e g scales, concluding that the presence of NLTE effects, 
and also a possible inconsistency between spectroscopic and 
photometric T c s scales, are probable explanations. Recent re- 
sults reporting offsets between spectroscopic and photomet- 
ric T c ff scales in cool stars, of similar magnitude, lend some 
credence to this interpretation. But we note that some au- 
thors claim consistency between the two scales, and that other 
sources of errors may be at play, such as uncertainties in the ab- 
solute flux calibration of photometric T^ffS- We also note that 
recent claims of such T c g offsets as caused by chromospheric 
activity cannot explain the present discrepancy given that both 
a Cen stars are considerably inactive stars. These discordant 
data still preclude a clear evaluation of the problem. For both 
a Cen A and B, the spectroscopic surface gravities agree well, 
within the uncertainties, with direct values derived from the 
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Fig. 10. The evolutionary state of the a Cen system. The stars 
are plotted in the HR diagram superimposed to the isochrones 
and evolutionary tracks of Kim et al. (2002) and Yi et al. 
(2003). The horizontal error bars refer to the uncertainties of 
Table[3] The actual uncertainties in the luminosities are smaller 
than the symbol size. The masses derived from evolutionary 
tracks, M A = 1.13 ± 0.01 and M B = 0.89 ± 0.05, agree very 
well with the orbital solution of Pourbaix et al. (2002). The 
tracks are labelled with masses in solar units, the numbers 
alongside the tracks are ages in Gyr. The thin solid lines be- 
tween the tracks join points with the same age. For aCen A, an 
age of 4.5 to 5.3 Gyr can be inferred. A match of aCen B within 
this age range is only possible if its T g is revised upwards (see 
text). 

dynamical masses of Pourbaix et al. (2002) and the radii of 
Kervella et al. (2003). The atmospheric parameters resulting 
from our analysis are collected in Tab [3] 

The abundance pattern of the system, when the various au- 
thors's data are considered for those elements for which at least 
three independent analyses are available, is found to be en- 
riched in Na, Mg, Si, Mn, Co, Ni and (with less reliability) 
Cu, and deficient in Y and Ba (Fig|9]l. This abundance pattern 
is found to be in very good agreement with recent results on 
the abundance ratios of metal-rich stars. Thus, the system may 
be considered as a normal pair of middle-aged, metal-rich, thin 
disk stars. 

An analysis of the evolutionary state of the system in the 
theoretical tracks of Kim et al. (2002) and Yi et al. (2003) 
yields a very good agreement of the evolutionary mass (Ma 
= 1.13 ± 0.01) and age (4.5-5.3 Gyr) of a Cen A with the re- 
sults of recent seismological and dynamical data (FigfTol). For 
a Cen B, a la upward revision of its T e g would bring its posi- 
tion in the HR diagram within reasonable agreement with the 
age found for a Cen A, and an evolutionary mass (Mb ~ 0.93) 
in good agreement with the dynamical one would result. This 
merely marginal compatibility suggests that, in order to fulfill 
the privileged situation of the a Cen system as a fundamental 
calibrator of the modelling of stellar structure and atmosphere 
models, additional analyses of component B seem to be neces- 



sary to quantify the onset and magnitude of possible NLTE in 
cool stars, as well as allow a more precise evaluation of possi- 
ble offsets between spectroscopic and photometric T e g scales 
in this class of objects. 
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